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ABSTRACT 

We simultaneously integrate in a numerical way the equations of motion of both the 
Magellanic Clouds (MCs) in the MOdified Newtonian Dynamics (MOND), MOdified 
Gravity (MOG) and Cold Dark Matter (CDM) frameworks for -1 < t < 1 Gyr 
in order to see if, at least in principle, it is possible to discriminate between them. 
Since the Large Magellanic Cloud (LMC) and the Small Magellanic Cloud (SMC) are 
at distances of approximately 50-60 kpc from the center of the Milky Way (MW), 
they are ideal candidates to investigate the deep MOND regime occurring when the 
characteristic MOND acceleration — 1-2 x 10~^° m s~^ is larger than the internal 
acceleration A of the system considered: indeed, the Newtonian baryonic accelerations 
involved are about 0.02 - O.OSAq for them. It turns out that CDM, MOND and 
MOG yield, in fact, different trajectories. In MOND also the External Field Effect 
(EFE) v4oxt must, in principle, be considered. Since for MW ^cxt ~ O.Ol^O) with a 
lingering uncertainty, we consider both the cases Acxt ^ ^Nj^cxt <C Aq and ^ext = 
AN,j4ext ^0- We also investigate the impact of the current uncertainties in the 
velocity components of MCs on their motions in the theories considered. In modeling 
the mutual interaction between the clouds and the dynamical friction (in CDM and 
MOND) we use for the masses of MCs the total (baryonic -I- dark matter) values, 
dynamically inferred, in CDM, and the smaller ones (baryonic), coming from direct 
detection of visible stars and neutral gas, in MOND and MOG. 
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1 INTRODUCTION 

In many astrophysical systems hke, e.g., spiral galaxies and 
clusters of galaxies a discrepancy between the observed kine- 
matics of some of their components and the predicted one on 
the basis of the Newtonian dynamics and the matter directly 
detected from the emitted electromagnetic radiation (visible 
stars a nd gas clouds) was present since the pioneeri ng stud- 
ies b: J ^ [Zwickv (I1933I) on the Coma cluster, and bv iBosmal 
jiJ) a 



(|l98ll) and lRubin et al.l (|l982l ) on spiral galaxies. More pre- 



cisely, such an effect shows up in th e galactic velocity rota- 
tion curves jPersic fc Salucc i 1996a b) whose typical pattern 
after a few kpc from the center differs from the Keplerian 
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^ He postulated the existence of undetected, baryonic matter; 
today, it is believed that the hidden mass is constituted by non- 
baryonic, weakly interacting particles in order to cope with cer- 
tain is sues pertaining galaxy formation and primordial nucleosyn- 
thesis llGondoldl2004l) . 



l/y'r fall-off expected from the usual Newtonian dynamics 
applied to the electromagnetically-observed matter. 

As a possible solution of this puzzle, the existence of 
non-baryonic, weakly-interacting Cold Dark (in the sense 
that its existence is indirectly inferred only from its grav- 
itational action, not from emitted electromagnetic radia- 
tion) Matter (CDM) was pro posed to reco ncile the pre- 
dictions with the observations l|Rubinl [19831 ) in the frame- 
work of the standard gravitational physics; for a general re - 
view on the CDM issue see, e.g, iKhalil fc Munoj l|2002h . 
while for the distr i bution of CDM in galaxies, see, e.g., 
ISalucci fc Borrielld (|2003l ). To be more definite, let us fo- 
cus on the Milky Way (MW) and adopt a very widely used 
model of its gravitational pot ential U . It consists of the stan- 
dard Miyamaoto-Nagai disk (jMivamoto fc Nagailll975l ) 

f^disk = ; ^ , (.1) 

the |Plumme3 l|l91ll ) bulge 



^bulgc 



CAfbulg. 

r + c 



(2) 
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and t he logarithmic CDM halo by iBinnev fc Tremaind 
l| 19871 ) 



+ d ) , (spherical halo) (3) 

with l|Law et aljlioosl : IWillett et al.ll2009l ) ^ = 1, k = 6.5 
kpc, b = 0.26 kpc, c = 0.7 kpc, Vhaio = 114 km s~^, 
d = 12 kpc. The m asses of the d isk and the bul g e use d 



by iLaw et al.l (|200d ) are those by IJohnston et"al] (|l999l ). 
i.e. Mdisk = 1 X 10" M0, Mbuigo = 3.4 x 10^° M© yield- 
ing a total baryonic mass of M = 1.34 x 10^ M©; how- 
ever, such a value is alm ost twice the mo st recent estimate 
(M = 6.5 x 10'° Mq) bv lMcGaughl | 20081 ) who includes the 

X 10'° Mp) and 



gas mass as well and yield Mdiak 
Mbuige = 2.07 X 10'° Mq. IXue et al.l l|2008l ) yield a total 
baryonic mass of M — 6.5 x 10'° M© as well; they use a 
different bulge-|-disk-|-CDM halo model of the Galaxy with 
Mdisk = 5 x 10'° Mr and Mbuige = 1.5 x 10^" Mq, as in 
ISmith et all (|2007l ). The model o f eg. Q-e q. (Et, with the 
laram eters' values bv lLaw et al.l ll2005h and I Johnston et all 
1999h . has been recentl y used bv IWillett et al.l ll2009l) to 
study the motion of the Icrillmair fc DionatosI ( 20061 ) tidal 



stellar stream a t Gala ctocentric distance of r < 16 — 18 kpc; 
I Read fc Moord (|2005l ) used it to study the motion of the 
tidal debris of the Sagittarius dwarf at 17.4 kpc from the 
center of MW. More specifically, the CDM halo model of eq. 
((3| corresponds to a CDM halo mass 



Mhalo 



G(r-2 + d?) ' 



(4) 



so that 



Mhaio(r = 60kpc)=3.5xl0'' Mq, (5) 

in agreement with the value bv lXue et al] (|2008| ) 

Mhaio('- = 60 kpc) = (4.0 ± 0.7) x 10^' Mq. (6) 

Concer ning Vhsio, othe r auth ors report different values for 
it; e.g., iRead fc Moord l|2005l ) use wo = 175 km s"\ where 
so that fhaio = 124 km s~' for them, while 



2w, 



IJohnston et al.l (|l999l ) yield the range 140 - 200 km s"' for 
their Ucirc = \/2«haio which maps into 70 km s~' < Uhaio < 
141 km s~'. However, it must be noted that values of Whaio 
too different from 114 km s~' would destroy the agreement 
of eq. (HI with the value of eq. (|6|l. 

Oppositely, it was postulated that the Newtonian laws 
of gravitation have to be modified on certain acceleration 
scales to correctly account for the observed anomalous kine- 
matics of such astrophysical systems without resorting to 
still undetected exotic forms of matter. One of the most 
phenomenologically successful modifications of the inverse- 
square Newtonian acceleration An, mainly with respect 
to spiral galaxies, is the Modified Newtonian Dynamics 
(MOND) (|Milgromlll983al lbll3) which postulates that for sys- 
tems experiencing total gra vitational acceleration A <^ Aq, 
with (iBegeman et al.lll99ll ) 



Ao 



(1.2 ±0.27) X 10"'° m s"^ 



^MOND = r. 



More precisely, it holds 



A 



m(X) 



, X: 



A 



(7) 



(8) 



(9) 



n{X) 1 for X ^ 1, i.e. for large accelerations (with re- 
spect to Ao), while /i(X) X yielding eq. ((8]) for X <^ 1, 
i.e. for small accelerations (again, with respect to Ao). The 
most widely used fo rms for the interpolating fun ction /i(X) 



most widely used to rms tor tne mterpolatmg tun ct 
are the "standard" (|Bekenstein fc Milgromlll984l ) 



X 



VT+x^' 

and the simpler (jPamaev fc Binne"vll2005h 

X 



i + x' 



(10) 



(11) 



It recently turned out that eq. (|11[) yields better results in fit- 
ting the terminal velocity curve of MW, t he rotation curve of 
the s t andard external galaxy NGC 3198 (iFamaev fc BinnevI 
l2005l : IZhao fc Famaevll20od : iFamaev et al.ll2007bl ) and of a 
sample of 17 high surface bright ness, early-type disc galaxies 
ISanders fc Noordermeerl [20071 ): eq. ^ becomes 



A = 



A^ 



l+\ 1 



4Ao 
An 



(12) 



with eq. ()lip. Eq. (|9]) strictly holds for co-planar, 
spherically and axially symmetric mass distributions 
iBrada fc Milgrom|[l99 5j); otherwise, the full modified (non- 
relativistic) Poisson equation (Bekenstoin & Milgrom 198^) 



V ■ 



\VU\ 
Ao 



VU 



4nGp 



(13) 



must be used. Attempts to yield a physical foundation 
to MOND, especially in terms of a relativistic covariant 
theory , can be be fo und i n, e.g.. iBekenstein fc Milgrom 



1984 ): iBekensteinI (|2004l ): iBruneton fc Esposito-Farese 
Mi 



20071 ): IZhao (2007^); for rece nt reviews of various aspect s 
of the MO ND pa r adigm, see Sanders fc McGaughl l)2002h : 



|Bek enstcin| (|2006h : iMilgromI l|2008l ). The compatibilitv of 
MOND with So la r System data has bee n investigated by 



MilgromI (Il983al): iTa lmadge et al.' (119881): ISereno _fc Jetzeil 



(I2006t ) : IBekenstein fc Maguciio (200^; ISanderd l|2006l ): 
lloriol (|2008al . |2009| ): IMilgromI (|2009i ). Generally speaking, 
many theoretical frameworks have been set up to yield 
a 1/r acceleration term able to explain the observed dy- 
namics of astrophysical systems; for exam ple, those en- 
compassing a logarithmic ext ra-potential l|Cadonil |2004| : 
iFabris fc Pereira Campo3l2009l) 



(7 = Cln I — 



(14) 



where C and r^, a length scale, are fit-for parameters. 
For other modified models of gravity used to explain, 
among other things, the galact ic rotation curve s with - 
out resorting to CDM see , e.g., ICapozziello et al . l|2006l) : 
iFrigerio-Martins fc SaluccH (|2007l ) : 'Moffat fc Toth' ('2008'). 

The Modified Gravity (MOG) (Moffat & Toth 2008) is 
a fully covariant theory of gravity which is based on the exis- 
tence of a massive vector field coupled universally to matter. 
The theory yields a Yukawa-like modification of gravity with 
three constants which, in the most general case, are run- 
ning; they are present in the theory's action as scalar fields 
which represent the gravitational constant, the vector field 
coupling constant, and the vector field ma ss. An approxi- 
mate solution of the MOG field equations l|Moffat fc TothI 
I2OO9I ) allows to compute their values as functions of the 
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source's mass. The resulting Yukawa-type modification of 
the inverse-square Newton's law in the gravitational field of 
a central mass M is 
G^M 



^MOG ~ 
Witlfl 



■{l + a[l-(l-^/ir)exp(-^r)]}r, (15) 



M 



[yM + E 
D 



Goo 

Gn 



M 

where Gn is the Newtonian gravitational constant and 



(16) 
(17) 



G 



E 
D 



20Gn, (18) 
25,000yMj^, (19) 
6,250^/M^kpc"\ (20) 

Such values have been obtained by iMoffat fc TothI l|2008l . 
I2OO9I ) as a result of the fit of the velocity rotation curves of 
some galaxies in the framework of the searches for an expla- 
nation of the rotation curves of galaxies without resorting to 
CDM. The validity of eg. I|15l) in the Solar System has been 
recently questioned in lloriol l|2008bl ). For (|McGaughl boOSl ) 
M = 6.5 X 10^° Mq, we have 



a ~ 16 
A = i: 



41 kpc. 



(21) 
(22) 



Traditionally, the phenomenology of both MOND and 
CDM paradigms is based on the electromagnetically de- 
tected matter (stars and gas cloud s) at no mor e than 
abo ut 20 kpc; i n view of the use bv IClewlev et al.l (|2004l ) 
and IXue et al] (|2008l ) of several recently discovered Blue 
Horizontal-Branch (BHB) stars as kinematical tracers at 
large radii (r ~ 60 — 130 kpc), it makes now sense to look 
at the remote periphery of the Galaxy as well to try to test 
CDM and alternative models of gravity. In this paper we 
wish to investigate the orbits of test particles at Galacto- 
centric distances r > 20 kpc, i.e. in the deep MONDian 
regime; we will use the Magellanic Clouds (MCs) moving at 
50-60 kpc from the center of MW. We will extend our anal- 
ysis also to MOG and to the action of CDM itself as well to 
see if our approach is able, at least in principle, to discrimi- 
nate between them; for another attempts on galactic scales, 
base d on the escape spee d in the solar neighborhood, see 
also iFamaev et al.l l|2007ar ) . At so large Galactocentric dis- 
tances many complications arising from an accurate mod- 
eling of the realistic distribution of mass can be avoided, 
both in MOND /MO G and in CDM frameworks. Mo r eover , 
iGardiner fc Noguchil (Il996l ). lYoshizawa fc Noguchil (|2003l ') 
and IConnors et al.l l|200 6^ demonstrated that the position 
of the Magellanic Stream (MS) follows the orbits of MCs. 
Therefore, it is interesting to compare the path of MS 
with the orbits predicted by CDM, MOND and MOG. 
Thus, it is hoped that our results will encourage more 
quantitative and detailed studies on MOND and MOG ap- 
plied to such systems; for numerical investigations on the 



2 iMoffat & TothI 1I2OO8I") used the eauivalent notation E c\ 



problem of the form ation of cosm ological str uctures and 
galact i c evolution, seelKnebe fc Gi bson ( 20 041) ; lHaghi et al.l 
12006"); Llinarcs ot al. ('2008D: lTiret fc Combeg lf2007l . l2008l ); 
Malekiani ct al.. (2009. '). 



2 MOTIONS IN CDM, MOND AND MOG: 
THE MAGELLANIC CLOUDS 

Concerning MOND and M OG, we will cons ider a central 
body with the same mass l|McGaughl [20081 ) M « 6.5 x 
IO^'^Mq of the total baryonic component of MW and a test 
particle distant several tens kpc from it, acted upon by the 
putative MOND/MOG gravitational fields of M. Such large 
distances allow to neglect the details of the real mass dis- 
tribution which may become relevant in MOND at closer 
distances (|Read fc Moorell2005l : iNipoti et al.ll2007h . To pre- 
liminarily test our approximation we applied eq. I|12p to 
the Sagittarius dwarf galaxy (r = Yl A kpc) and confronted 
the numerically integrated orbital sections i n the coordinate 
planes of the trajectory to those obtained bv lRead fc Moor3 
(2005') by using a no n-pointl ike baryonic potential (upper 
panel of Figure 2 in (|Read fc Moore 2005,)); we used the 
same integration interval of — 1 Gyr < t < 1 Gyr and the 
same baryonic mass (M = 1.2 X 1O"M0) bv lRead fc Moord 
(200^. It turns out th at we were succ e ssful in reproducing 
the orbital sections by iRead fc Moor3 Hooi); thus, we are 
confident of the validity of our approximation for the larger 
Galactocentric distances we will use in the following analy- 
sis. 

Another issue which, in principle, should be taken into 
account in MOND is the so-called External-Field-Effect 
(EFE); it may become relevant wit h cluster of galaxies 
Wu et al.l 120071'). According t o. e.g.. ISanders fc McGaughl 



20021 ) ; lFamaev et al.l (|2007al ): [Angus fc McGaughl (|2008l ) 



+ A\ 



Ao 



(23) 



where An is the Newtonian acceleration of the system alone, 
A is its total internal acceleration, while Aoxt denotes the ac- 
celeration induced by any external field. By using the sim- 
pler form of eq. (|ll|l for ^, one approximately obtains from 
eq. nisi 



A: 



An 



1 - 



Ak 



4Ao / ^ 
An V Ao 



(24) 

For Ao 0, A ^ An, as expected. For A^xt — > 0, i.e. 
Aaxt «C Ao and Acxt «C An, one has A eq. (|12|l . For 



Acxi 



< 1 



only, the total acceleration becomes 
An 



A', 



while for 



1 



1 - 



+ 



4Ao 
An 



(25) 



(26) 



A 
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Table 1. Larg e Magella nic Cloud (LMC): coordinates 
llKallivavalil et al. J2OO6; W u et alj boOSh . in kpc, and veloc- 
itv components l|Wu et al.l boOSh . in km s~l, of LMC in a 
Galactocentric rest frame {X, Y, Z} with the Z— axis point- 
ing toward the Galactic north pole, the X— axis pointing in 
the direction from the Sun to the Galactic center, and the 
y— axis pointing in the direction of the Suns Galactic rotation 
llKallivavalil et a,l.ll2006l : iBesla et al.ll2007l) . They yield r = 49.5 
kpc, V = 378 km s~^. The un certainties in the coordinates can 
be neglected llCioni et al.ll200d h. 



Xo 
Xo 



-0.8 



i ± 12 



Yo 
Yo 



-41.5 
-268 ± 11 



Zo = -26.9 
Zo = 252 ± 16 



Table 2. 



Small Magella nic 
[2006; 



C loud (SMC): coordinates 
llKallivavalil et al.|j 2006: Wu et aD I2OO8I ). in kpc, and veloc- 
ity components 11 Wu et al.. ,20081) . in km s~^, of SMC in a 
Galactocentric rest frame {X, Y, Z} with the Z— axis point- 
ing toward the Galactic north pole, the X— axis pointing in 
the direction from the Sun to the Galactic center, and the 
y — a xis point ing in the direction of the Suns Galactic rotation 
llKall ivavalil ct al.ll2006l : iBesla et allbOOTl ). They yield r = 58.9 
kpc, V = 301 km s~^. The uncertainties in the coordinates are 
negligible l|Cioni et al.ll2000h . 



Xo = 15.3 

Xo = -87 ± 48 



Yo 
Yo 



-36.9 
-247 ±42 



Zo = -43.3 
Zo = 149 ± 37 



only, it is 



A: 



AnAo 1 + 



Interestingly, if 



An 



'a'o 



« 1, 



then 



A ^ VAj^Ao = 



VGMAo 



(28) 



(29) 



(30) 



In the case of MW, it is very difficult to reliably assess the ex- 
ternal field because it may be due to several factors like, e.g., 
the Large Scale Structure and the Great Attractor region 
(Acxt/Ao = 0.01), but also the galaxy M31 Andromeda, at 
800 kpc from MW, and the Coma and Virgo clusters, whose 
field are t i me- va rying, may play a role. For a discussion see 

et al.l (|2008l ) . In view of the lingering uncertainty of Acxt , 
in the following we will use eq. I|12|l : however, we will also in- 
vestigate the case in which >lext = Ajsi, ^lext ^ because 
it may occur in MW at the large Galactocentric distances 
considered here. 

As a concrete example of motion in deep MOND regime 
(4nMo ^ 0.03 - 0.02), let us consider both MCs; the Large 
Magellanic Cloud (LMC) is at 49.4 kpc from the center 
of MW (Galactic Center, GC), while the Small Magellanic 
Cloud (SMC) is located at 59 kpc from GC. LMC and 
SMC's Galactocen tric cartesian coordinates and velocities 
l|KallivavaUl et al.1 2006: Wu et al...2008l ) are in Table [1] and 
Table [S] It can be noted that the velocity components of 
LMC are uncertain at more than 4 — 14%. The situation 
for the position components is much better since they are 
known with uncertainties in the range 0.1 — 1%, as it results 



from the analysis o f the tip of the red giant branch (TRGB) 
applied to MCs bv lCioni et al] l|2000l ): thus, we will neglect 
them in the following. Also for S MC the un certainty in the 
position components is negligible (|Cioni et al. 2000). while 
the velocity components are known at 10%. 

We simultaneously integrated in a numerical way the 
equations of motion of both MCs in MOND, MOG and CDM 
by using the initial conditions of Table [T] and of Table [2] for 
— 1 < t < 1 Gyr. In addition to the main pull due to MW, we 
also included the mutual attractions of MCs and the effect 
of the dynamical fri ction due to their motion through the 
Galactic dark halo (|Binnev fc Tremainelll987|): the mutua l 
dynamical friction was neglected ( Kallivavalil et al]|2006l ). 
Concerning the pull by LMC on SM C, we mo d eled i ts ac- 
tion in N ewtonian dynamics from a iPlummerl |l9l3)-type 
potential (|Kallivavalil et al.|[2006l ) 



f^LMC = 



Gm 



LMC 



2LMC)^ + -K'lmC ' 



(31) 

with i^LMC ~ 3 kpc. In MOND, since the accelera- 
tion imparted by LMC on SMC is of the order of about 
0.05j4o, we adopted eq. (jSj with M ttilmc and 

r = a;LMc)^ + {y- VLMc)^ + {z- 2LMC)^ + -K'EmC' 

while in MOG we used eq. (|15|) with Ad — * mLMC and 
r = a^LMc)^ + (y- VhMc)^ + {z- zlmc)^ + -K'Emc- 

An analogous expression fo r the pull by SMC on LM C holds; 
in this case, A'sMC = 2 kpc (|Kallivavalil et al.l200d ). The dy- 
namical friction experienced by, say, SMC in going through 
the dark halo of the Galaxy has been modelled, in CDM, as 



D 



V 

^fric 



with (|Kallivavalil et aLllioO^ : iHaghi et al.ll2009l ) 



0.428 In A 



Gmsh 



(32) 



(33) 



where the Coulomb logarithm In A « 3 (|Binnev fc Tremainel 
Il987h . We also included the mutual dynamical friction ex- 
perienced by SMC when its distance from LMC gets smaller 
than 15 kpc jBekki fc Chiball2005l : iKallivavalil et al.ll2006l ) 
by replacing in eq. (|33p InA = 3 with IuAls ~ 0.2 and r 
with Tmutuai- The dynamical friction plays a non-negligible 
role also in several astrophysical systems in the framework of 
MOND dCiotti fc Binnevll2004l : lsrnchez -Salcedo et al. I I2OO6I : 



iNipoti et al.l 



( Ciotti fc Binnevlf2004l : iNipoti et al.ll2008h 



2008f): in our case, we mod el it by assuming 



(34) 



• MOND 
'^fric 

''fric 



V2 



V2 



1 + 



A 
An 



Since a model of the dynamical friction has not yet been 
developed in the framework of MOG, we did not include it. 

Concerning MCs' masses entering both their mutual in- 
teractions and the dynamical friction, for consistency rea- 
sons we adopted the t otal (baryonic -I- dark matter) values 
mLMC = 2 X 10^°Mpi ([ Kallivavalil et al." 2006 ; HaEhi et alj 
i2009i : ISchommer et al.ni992 ; Gardiner fc Noguchi 1996^ 
coming from radial velocities of several of the oldest star 
clusters in LMC lying well beyo nd 6 kpc from its cen- 
ter, andfl msMC = 3 x IO^Mq (jKallivavalil et all I2OO6I : 

^ This typical value has been chosen by IKallivavalil et all 
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iHaghi et all |2009| ) when integrating the CDM model. In- 
stead, we used the smaller, baryonic values mLM C = (2.7 + 
0.5 = 3 .2) X lO^Mp) (vi sible disk + neu tral gas ()Kim et all 
1 19981 )) (|van der Marel et al. 20oj . l2008l ) and msMC = (3.1 + 
5.6 = 8.7) X IO^Mq (to tal stellar mass + neutral gas) 
(|van der Marel et al.ll2008l ) for MOND and MOG. For a re- 
cent discussion of the method s employed to obtain such fig- 
ures and of other results, see l|van der Marel et al.1120081 ). 

2.1 The Large Magellanic Cloud 

In Figure [T] we show the sections in the coordinate planes 
of the LMC's orbits for CDM (red dash-dotted curves), 
MOND (blue dashed lines and light blue dotted lines) and 
MOG (yellow continuous curves) over — 1 < t < 1 Gyr. 
The dynamical models and their parameters' values are 
those described in Section [l We used eq. for MOND 
(^ext <K and ^ext <C ^o) obtaining the blue dashed lines 
depicted; indeed, for LMC A^/Ao = 0.03, so that eq. (fT2)) is 
adequate for it by assuming ylext = O.OIAq. Concerning the 
impact of EFE in MOND on LMC, we also investigated it 
in the case Aext = and Acxt Aq; thus, we numerically 
integrated trajectories with eq. (|30|) as well, which corre- 
sponds to an external field equal to the internal Newtonian 
one, obtaining the light blue dotted curves shown. The same 
approach will be used in Section [2.21 for SMC. The middle 
panel of Figure [5] shows the Galactocentric distance of LMC 
for the central values of the velocity components of Table [1] 
over —l<t< 1 Gyr. 

The smallest Galactocentric distance occurs for MOND, 
while MOG and CDM yield the largest one amounting to 
about 255-270 kpc after +1 Gyr and 250-280 kpc after -1 
Gyr. MOG and CDM differ by about 15 kpc, while the dis- 
crepancy between MOG/CDM and MOND is approximately 
70 — 80 kpc after —1 Gyr; in the past Gyr the discrepancy 
between MOG and CDM is of the order of 30 kpc, while 
MOND differs from MOG/CDM by a few 10 kpc. Over the 
next Gyr the Galactocentric distance of LMC undergoes a 
steady increase. I t may be interesting to recall that, accord- 
ing to IWu et al] (|2008l ). LMC is on a bound orbit; how- 
ever, they did include neither the mutual interaction with 
SMC nor the dynamical friction. The difference between the 
MOND trajectories for ^ = X/ (1 -I- X) and ^ = X \s rather 
small; discrepancies of the order of 10 kpc or less occur at ±1 
Gyr. All the models considered tend to undergo reciprocal 
departures after some about ±500 Myr. 

Figure [3] shows the impact of the dynamical friction; 
after +1 Gyr, without modeling it in CDM and MOND, the 
mutual difference between CDM and MOG tends to increase 
by about 10 kpc, while the MONDian trajectories are left 
almost unaffected. Instead, at —1 Gyr the discrepancy be- 
tween CDM and MOG gets reduced by 10 kpc, while the 
MOND distance is smaller by about 20 kpc. 

The impact of the uncertainties in the velocity compo- 
nents of LMC has been evaluated as it will be done for SMC 
in Section 12.21 it is shown in the upper (minimum velocity) 

after an examination of the values coming from obser- 
vations of carb on stars d Hardv et al.|[l989l ) and planetary nebulae 
llDopita et al.l llQSS). and from a virial analysis of the kinematic s 
of thousands of red giant stars in SMC llHarris fc Zaritskvll2006f) . 
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Figure 1. Sections in the coordinate planes of the numerically 
integrated trajectories of LMC experiencing: a) The Newtonian 
acceleration with CDM (red dash-dotted line) b) The MOND ac- 
celeration with ^ = X/(l + X) (blue dashed line) c) The MOND 
acceleration with fi = X (light blue dotted line) d) The MOG 
acceleration (yellow continuous line). The central values of the 
initial conditions of Table [T] have been used. For th e baryonic 
masse s of MW's bulge and disk we used the values bv lMcGaughl 
||2008|) . with a total baryonic mass oi M = 6.5 X 10^° Mq. For the 
masses of MCs, entering their mutual interactions and the dynam- 
ical friction, both modelled in the present integration, the total 
values (baryonic -|- dark matter) have been adopted for CDM, 
while those encompassing only baryonic components have been 
used for MOG and MOND. The time span of the integration is 
-1 < t < 1 Gyr. 
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Figure 2. LMC: Galactocentric distance r, in kpc, for —1 < 
* < 1 Gyr. Red dasli-dotted line: CDM. Blue daslied line: MOND 
(/i = XI{\^X)). Liglit blue dotted line: MOND (^t = X). Yellow 
continuous line: MOG. The initial condition for the position is 
r = 49.4 kpc. Upper panel: for the velocity we adopted = 
-86 -I- 12 = -74 km s-\ yo = -268 + 11 = -257 km s"!, 
20 = 252 — 16 = 236 km s" ^ yielding the minimum value v = 
356.7 km s~^. Middle panel: the central values of Table [T] have 
been adopted for the velocity. Lower panel: for the velocity wc 
adopted io = -86 - 12 = -98 km s"!, yo = -268 - 11 = -279 
km s" ^, io = 252 + 16 = 268 km s" ^ yielding the maximum 
value V = 399.1 km s~^. For the masses of MCs, entering their 
mutual interactions and the dynamical friction, both modelled in 
this integration, the total values (baryonic + dark matter) have 
been adopted for CDM, while those encompassing only baryonic 
components have been used for MOG and MOND. 



and lower (maximum velocity) panels of Figure |2l Differ- 
ences with respect to the nominal case are present. Indeed, 
for the smallest value of the velocity (io = —86 + 12 = —74 
kms"\ yo = -268-Hl = -257kms~\ io = 252-16 = 236 
km s~^ ), the overall discrepancy among CDM/MOG and 
MOND is of the order of 50 - 60 kpc after -1-1 Gyr, with a 
reduction of the final distances in CDM/MOG with respect 
to the middle panel of Figure [2] (20 kpc for CDM, 30 kpc 
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Figure 3. LMC: Galactocentric distance r, in kpc, for —1 < t < 
1 Gyr. Red dash-dotted line: CDM. Blue dashed line: MOND 
= X/{l + X)). Light blue dotted line: MOND {fi = X). Yellow 
continuous line: MOG. The initial condition is r = 49.4 kpc for the 
position; for the velocity, the central values of Table [T]have been 
adopted. No dynamical friction has been applied in CDM and 
MOND. For the mass of SMC the total value (baryonic + dark 
matter) has been adopted for CDM, while that encompassing only 
baryonic components has been used for MOG and MOND. The 
time span of the integration is — 1 < t < 1 Gyr. 



for MOG and 50 kpc for MOND); after -1 Gyr the CDM 
distance is 30 kpc smaller than for the nominal values of the 
velocity components, MOG is 35 kpc below the level of the 
middle panel of Figure [2j while the MOND curves experi- 
ence a reduction of about 40 kpc. For io = —86 — 12 = —98 
kms"\ 2/0 = -268-11 = -279 kms"\ io = 252 + 16 = 268 
km s~^, corresponding to the maximum velocity, the rela- 
tive discrepancy after -1-1 Gyr among the various models is 
about 60 — 80 kpc, with an increase of each of them with 
respect to the middle panel of Figure [2] (25 kpc for CDM, 
30 kpc for MOG, and 50 - 60 kpc for MOND); also at -1 
Gyr there is an overall increase with respect to the case of 
the nominal values of the velocity components (20 kpc for 
CDM, 30 kpc for MOG, and 30 - 40 kpc for MOND). 



2.2 The Small Magellanic Cloud 

Figure |4] depicts the sections in the coordinate planes of the 
SMC's orbits over -1 < t < 1 Gyr for CDM (red dash- 
dotted curves), MOND (blue dashed lines and light blue 
dotted lines) and MOG (yellow continuous curves) by using 
the central values of the initial velocities of Table [2] and the 
same values of Section 12.11 for the masses of MW and MCs 
and of the other models' parameters. Concerning MOND 
and the impact of EFE, we followed the same approach as 
for LMC in Section [2T] It can be noted that MOND, MOG 
and CDM yield different orbital patterns, especially in the 
{xy} and {xz} planes, and it is possible, in principle, to 
discriminate among them. 

In the middle panel of Figure [5] we plot the time evo- 
lution of the Galactocentric distance of SMC according to 
CDM, MOND and MOG in the next Gyr for the central 
values of the velocity components of Table (2] The distance 
reached in all the three models after -1-1 Gyr is practically 
the same, amounting to 220 — 230 kpc. After —1 Gyr the 
scatter among the models considered is larger, amounting to 
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Figure 4. Sections in the coordinate planes of the numerically 
integrated trajectories of SMC experiencing: a) The Newtonian 
acceleration with CDM (rod dash-dotted line) b) The MOND 
acceleration with ^l = X/(l+X) (blue dashed line) c) The MOND 
acceleration with fi = X (light blue dotted line) d) The MOG 
acceleration (yellow continuous line). The central values of the 
initial conditions of Table |2] have been used. For th e baryonic 
masse s of M W's bulge and disk we used the values by iMcGaughl 
||2008| '). with a total baryonic mass of M = 6.5 X lO^" Mq. For 
the masses of MCs, entering their mutual interactions and the 
dynamical friction, both modelled in this integration, the total 
values (baryonic -|- dark matter) have been adopted for CDM, 
while those encompassing only baryonic components have been 
used for MOG and MOND. have been adopted. The time span of 
the integration is — 1 < t < 1 Gyr. 
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Figure 5. SMC: Galactocentric distance r, in kpc, for —1 < 
t < 1 Gyr. Red dash-dotted line: CDM. Blue dashed line: MOND 
i^i = X/{l + X)). Light blue dotted line: MOND (^t = X). Yellow 
continuous line: MOG. The initial condition for the position is 
r = 58.9 kpc. Upper panel: for the velocity we adopted xq = 
-87 -I- 48 = -39 km s"!, yo = -247 + 42 = -205 km s"!, 
io = 149 — 37 = 112 km s" yielding the minimum value v = 236 
km s~^. Middle panel: the central values of Table |2] have been 
adopted for the velocity. Lower panel: for the velocity we adopted 
±0 = -87-48 = -135 km s"!, yo = -247-42 = -289 km s"!, 
io = 149-1-37 = 186 km s" yielding the maximum value v = 369 
km . For the masses of MCs, entering their mutual interactions 
and the dynamical friction, both modelled in this integration, 
the total values (baryonic + dark matter) have been adopted for 
CDM, while those encompassing only baryonic components have 
been used for MOG and MOND. 



about 40 kpc. As for LMC, all the models considered tend to 
undergo reciprocal departures after some about ±500 Myr, 
although smaller. 

Figure [6] shows that switching off the dynamical friction 
in CDM and MOND does not substantially alter the overall 
picture. 

The uncertainty in the velocity components of SMC 
may have different consequences on its orbit for the mod- 
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Figure 6. SMC: Galactocentric distance r, in kpc, for —1 < t < 
1 Gyr. Red dash-dotted line: CDM. Blue dashed line: MOND 
(fi = X/{l + X)). Light blue dotted line: MOND (^t = X). Yellow 
continuous line: MOG. The initial condition is r = 49.4 kpc for 
the position; for the velocity, the central values of Table |2]have 
been adopted. No dynamical friction has been applied in CDM 
and MOND. For the mass of LMC the total value (baryonic + 
dark matter) has been adopted for CDM, while that encompassing 
only baryonic components has been used for MOG and MOND. 



els considered. The Galactocentric distance of SMC for the 
maximum value of its speed, i.e. wsmc = 369 km s~^ corre- 
sponding to io = -87-48 = -135 kms~\?/o = -247-42 = 
-289 km s~\ io = 149 -f- 37 = 186 km s"\ is shown in the 
lower panel of Figure [5] By comparing it with the middle 
panel of Figure [5] it can be noted that, after -1-1 Gyr, the 
Galactocentric distance increases, in particular in CDM and 
MOG; instead, at —1 Gyr the increase is more uniform for all 
the models. In the upper panel of Figure[5]we depict the case 
for ±0 = -87-^48 = -39 km s-\ yo = -247+42 = -205 km 
s~^, io = 149 — 37 — 112 km s~^ yielding the minimum value 
for the SMC's speed v = 236 km s~^. In this case, CDM and 
MOG yield a smaller Galactocentric distance after -|-1 Gyr: 
indeed, it is as large as 140— 150 kpc. Instead, after —1 Gyr 
the MOND curves lie in between the MOG/CDM ones, with 
an overall reduction of 50 — 60 kpc for all the models. 



2.3 The mutual distance between SMC and LMC 

In the middle panel of Figure [7] the 

m utual SMC-LMC distance A = 

\/ (a;LMC — a;sMc)^ + (j/lmc — ysMc)^ + (zlmc — zsmc)^ is 
shown for the central values of the velocity components of 
both MCs an with the same models and parameters' values 
of Section (2] The pattern by CDM is quite different with 
respect to those by MOND and, to a lesser extent, MOG, 
both in the size of the distance reached and, especially, in 
the temporal signature. After -1-1 Gyr, MOND exhibits a 
bounce yielding the smallest maximum reciprocal separa- 
tion, i.e. about 25 kpc, while for CDM, which yields an 
increasing signal, it is approximately 50 kpc. Instead, at —1 
Gyr CDM and MOG reach 80 kpc, while MOND is around 
20 kpc. 

In Figure [S] we show MCs mutual distance without dy- 
namical friction. After -f-l Gyr, the CDM maximum distance 
is 60 kpc, while the MOND curves tend to approach the 



MOG one at 30 kpc. After -1 Gyr, CDM reaches about 70 
kpc; MOND is below the 20 kpc level. 

The impact of the uncertainties in the velocity of both 
SMC and LMC on the mutual separation is depicted in the 
upper and lower panels of Figure [T] Low velocities (upper 
panel) yield a change of the bouncing time for MOND and 
an increase of A in CDM and MOG by more than 50 kpc 
for ±1 Gyr. Note also the increase of the MOND curves at 
— 1 Gyr: they pass from 20 kpc to 60 kpc, while they are not 
substantially changed after -1-1 Gyr. On the contrary, high 
velocities (lower panel) tend to yield an overall reduction 
of the distances among all the models, with Acdm reduced 
down to 40 kpc. 

Generally speaking, the repeated close encounters (in 
MOND and MOG) may have an impact on the star for- 
mation history and/or morphology of both MCs; anyway, 
discussing such interesting issues is beyond the scope of the 
present paper. 



3 SUMMARY AND CONCLUSIONS 

We simultaneously integrated in a numerical way the orbits 
of both MCs for -1 < t < 1 Gyr within MOND, MOG 
and CDM to see if, at least in principle, it is possible to 
discriminate among them. This is, in principle, important 
also because it is believed that MS follows the orbits of MCs. 

Since LMC and SMC are at about 50-60 kpc from the 
Galactic center, they are ideal candidates to investigate the 
deep MOND regime {Aisi/Ao = 0.03 - 0.02); moreover, the 
details of the realistic mass distribution can be neglected. 
Thus, for MOND (and MOG) we used a pointlike approxi- 
mation for the baryonic mass of MW; we tested it by success- 
fully reproducing the orbital paths of the Sagittarius dwarf 
galaxy (r = 17 kpc) obtained by other researchers with the 
MONDian fully non-linear modified Poisson equation. For 
CDM we used a logarithmic halo potential which is able to 
reproduce the value of the Galactic mass at 60 kpc obtained 
independently by analyzing different tracers. We also took 
into account the mutual MCs interaction and, for MOND 
and CDM, the dynamical friction as well. For the masses 
of MCs we used the total (baryonic -I- dark matter) values 
dynamically inferred in CDM, and the smaller ones com- 
ing from direct detection of the electromagnetic radiation 
emitted by stars and neutral gas in MOND and MOG. 

It turns out that, in fact, CDM, MOND and MOG do 
yield different trajectories for SMC and LMC. In general, 
the spatial extension of the orbits' sections in the coordi- 
nate planes is larger for CDM and MOG with respect to 
MOND. SMC experiences larger discrepancies among the 
various models than LMC. Since for MW ^cxt ~ O.OlAo, 
we also investigated EFE in MOND by considering not only 

Aoxt < AN,Aoxt < Ao, but also Aext ~ AN,Acxt < Aq 

which cannot be excluded in view of the lingering uncer- 
tainty in MW's EFE. The resulting orbital patterns are 
rather similar to those obtained neglecting Acxt- We also 
investigated the impact of the present-day uncertainties in 
the velocity components of MCs on their trajectories in the 
models considered by finding that it is more notable for SMC 
than LMC; CDM and MOG are more sensitive to such a 
source of bias than MOND. In general, the largest discrep- 
ancies among the various models occur around ±1 Gyr. This 
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suggests that extending the integration time may yield in- 
teresting findings; it may be the subject of further analyses. 
Over the timescale considered, the dynamical friction does 
not make the paths too much different in the various models 
examined. 
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Figure 7. Mutual distance A between SMC and LMC, in kpc, for 
-1 < t < 1 Gyr. Red dash-dotted line: CDM. Blue dashed line: 
MOND (^t = Light blue dotted line: MOND (/.t = X). 

Yellow continuous line: MOG. Upper panel: for the velocity com- 
ponents of both SMC and LMC we used their minimum values. 
Middle panel: for the velocity components of both SMC and LMC 
we used their central values. Lower panel: for the velocity compo- 
nents of both SMC and LMC we used their maximum values. For 
the masses of MCs, entering their mutual interactions and the dy- 
namical friction, the total values (baryonic + dark matter) have 
been adopted for CDM, while those encompassing only baryonic 
components have been used for MOG and MOND. 
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Figure 8. Mutual distance A between SMC and LMC, in kpc, 
for -1 < t < 1 Gyr. Red dasli-dotted line: CDM. Blue dashed 
line: MOND (^t = X/(l + X)). Light blue dotted line: MOND 
{/i = X). Yellow continuous line: MOG. For the velocity compo- 
nents of both SMC and LMC we used their central values. No 
dynamical friction has been modeled. For the masses of MCs, en- 
tering their mutual interaction, the total values (baryonic -|- dark 
matter) have been adopted for CDM, while those encompassing 
only baryonic components have been used for MOG and MOND. 



